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Abstract. We report the abundance determination in the atmosphere of the bright halo star HD221170. The 
spectra were taken with the Terskol Observatory's 2.0-m telescope with a resolution R=45000 and signal-to-noise 
ratio up to 250 in the wavelength region 3638-10275 A. The adopted atmospheric parameters correspond to an 
effective temperature Tcff=4475 K, a surface gravity log g—1.0, a microturbulent velocity Umicro=1.7 km s~^, and 
a macroturbulent velocity Umacro=4 km s~^. The abundances of 43 chemical elements were determined with the 
method of spectrum synthesis. The large overabundances (by 1 dex relative to iron) of elements with Z> 38 are 
shown to follow the same pattern as the solar r-abundances. The present HD221170 analysis confirms the non- 
universality of the r-process, or more exactly the observation that the astrophysical sites hosting the r-process do 
not always lead to a unique relative abundance distribution for the bulk Ba to Hg elements, the Pb-peak elements, 
and the actinides. 
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1. Introduction 

The chemical composition of different objects, particularly 
halo stars, is one of the most important clues for un- 
derstanding the structure and evolution of the Universe. 
Typical investigations of the chemical compositions of 
stars usually show 15-40 chemical elements, depending 
the type of the star. The most detailed stellar abundance 
patterns consist of 50-58 elements (see Yushchenko et al. 
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2004 for examples). The best example is the abundance 
pattern of the halo star CS22892-52 with the determina- 
tion of 58 elements (Sneden et al. 2003). This star and 
another three r-process-rich stars CS31082-001 (Cayrel et 
al. 2001, HiU et al. 2002 ), HD115444 (Westin et al. 2000) 
and BD-f 17°3248 (Cowan et al. 2002) are halo stars with 
known enhanced abundances of thorium with respect to 
iron. In the specific case of CS31082-001, the determina- 
tion of the uranium abundance also enabled investigators 
to derive from the thorium to uranium abundance ratio 
a stellar age of 13±4 billion years. Such an estimate di- 
rectly provides an independent lower limit for the age of 
our Galaxy and the Universe. 
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The investigation of the thorium abundance in the at- 
mosphere of different stars of our Galaxy has a long his- 
tory. Buchcr (1987), Morel et al. (1992) and Francois et 
al. (1993) studied a set of disk and halo stars using the 
strongest Th line at 4019.129 A. Yushchenko & Gopka 
(1994) determined the thorium abundance in Procyon 
using 4 faint lines in the 3200-3500 A spectral range, 
while Gopka et al. (1999) determined Th at the surface of 
Arcturus using other lines. 

A special effort was later dedicated to the observation 
of an increasing number of Th-rich halo stars (see for ex- 
ample Johnson & Boltc, 2001). But the number of lines 
considered for the abundance determination did not nec- 
essarily increase; for example the latest investigation of 
Honda et al. (2004) is based on the single 4019 A Hne. 
Numerous thorium lines were detected only in the four 
above-mentioned Th-rich stars. 

Many authors have used the Th/Eu and other ratios to 
estimate the age of the observed stars, although Goriely & 
Clerbaux (1999) and Goriely & Arnould argued that the 
Th-to-Eu abundance ratio is not a reliable chronometer to 
derive the stellar age. New observations by Honda et al. 
(2004), as well as the present analysis confirm this result. 

HD221170 is one of the bright template halo stars. 
This seven-magnitude star has been investigated regularly 
as new methods and new observational facilities became 
available. The first paper in this series was Wallerstein et 
al. (1963), followed by Gilroy (1988) and many others. The 
last detailed abundance pattern of heavy elements in this 
star was provided by Burris et al. (2000) . Yushchenko et al. 
(2002) published preliminary results for the abundances of 
elements heavier than Dy. However, for such elements the 
majority of the spectral lines are located at wavelengths 
shorter than 4500 A, and the low signal-to-noisc ratio ob- 
tained in this spectral region at that time required fur- 
ther investigations. Gopka et al. (2004; hereafter Paper I) 
used the spectral data from the above-mentioned work 
and the spectrum from the archive of the Hautc-Provence 
Observatory 1.88 m telescope to derive the abundance of 
elements with Z < 68. Previous investigations concerning 
the surface abundances of HD221170 and the determina- 
tion of atmospheric parameters are reviewed in Paper I. 
According to Burris et al. (2000), Gopka et al. (2001) and 
Paper I, the iron abundance in HD221170 is deficient by 
about 2 dex with respect to the sun while the abundance 
of elements heavier than barium are overabundant with 
respect to iron. Nevertheless, since the spectrum used in 
Paper I from the Haute-Provence Observatory covers the 
limited wavelength range of 4480-6820 A, we were unable 
to determine the abundance of elements heavier than Dy. 

In Sect. 2, new observations of higher quality and 
broader wavelength coverage are detailed and the adopted 
atmospheric parameters described. In Sect. 3 the result of 
the abundance analysis is presented. In Sect. 4, the sur- 
face abundances of HD221170 are compared with those of 
the other r-process-enriched halo stars already observed, 
and their implication concerning the non-universality of 
the r-process is discussed. 



2. Observations and atmospheric parameters 

The observations were conducted with the Coude-echelle 
spectrometer (Musaev et al. 1999) mounted at the 2-m 
" Zeiss" telescope at the Peak Terskol observatory located 
near Mt. Elbrus (Northern Caucasus, Russia) 3,124 m 
above sea level. The spectrometer is used in the mode 
with a spectral resolution of 45,000. The best signal-to- 
noise ratio is over 250 in the red spectral region and the 
wavelength coverage includes the 3638-10275 A range. 
Galazutdinov (1992) DECH20 and Yushchenko (1998) 
URAN codes are used for the data processing. 

In Paper I, based on previous investigations of the 
star and new spectral data, the following model atmo- 
sphere parameters of HD221170 were found to be op- 
timum: Teff=4475±50 K for the effective temperature, 
log g=1.0±0.1 for the gravity, Wniicro=l-7±0.1 km s"'^ for 
the microturbulence and t;macro=4-0±0.5 km s~^ for the 
macroturbulence. 

These parameters and the abundances from Paper I 
are used to construct individual atmosphere models using 
the standard Kurucz's ATLAS12 code. We find for the 
final model that the correlation(s) of iron abundances with 
excitation potentials and equivalent widths of iron lines 
are close to zero. 

We also calculated the models corresponding to the pa- 
rameter sets (reff=4475 K; log 5=1.2) and (Teff=4575 K; 
log gf=1.0). These models are used to estimate the abun- 
dance errors due to possible uncertainties in the adopted 
value of the temperature and gravity. 

3. Abundance analysis 

The equivalent widths of the iron lines are estimated 
by fitting their profiles with a Gaussian function. The 
iron abundance is calculated using the model atmosphere 
method on the basis of the Kurucz (1995) WIDTH9 pro- 
gram. In contrast, (iiffcrc;ntial spcc;trum synthesis meth- 
ods are used for all the other elements. For each line, we 
tried to find its counterpart in the solar spectrum atlas of 
Delbouille et al. (1973). This procedure frees us from un- 
certainties connected with oscillator strengths of spectral 
lines. The URAN code (Yushchenko 1998) and SYNTHE 
spectrum synthesis program (Kurucz 1995) are used to 
approximate the observed spectrum with the synthetic 
one. The solar abundances of Grevesse & Sauval (1999) 
arc considered. A synthetic spectrum of HD221170 for 
the whole wavelength range helps us to identify spec- 
tral lines. SYNTHE program (Kurucz 1995) was used to 
produce the synthetic spectrum. Atomic and molecular 
lines are included from Kurucz (1995) as well as Morton 
(2000), Biemont et al. (2002) and partially from the VALD 
database (Piskunov et al. 1995) 

Holweger's partition function for thorium is used 
(Morell et al. 1992). Hyperfine structure and isotopic split- 
ting are taken into account for Li, Sc, V, Mn, Co, Cu, Ba, 
and Eu. The splitting data for Li are taken from Shavrina 
et al. (2003), for Ba - from Francois (1996) and for other 




Fig. 1. The observed spectrum of HD221170 (squares) 
and the synthetic spectra (soUd Unes) calculated with 
our final abundances. The axes are the wavelength in 
angstroms and relative fluxes. The positions of the spec- 
tral lines taken into account in the calculations are marked 
in the bottom part of the figure by short and long dashes. 
For some of the strong lines, the identification and isotopic 
composition for molecular lines are given. The position of 
the Hf II 4093.155 A line is marked by a vertical dotted 
line. The different synthetic spectra correspond to a Hf 
abundance lower or higher by 0.5 dex with respect to the 
abundance obtained from the optimum value. 

elements from Kurucz (1995). We found counterparts in 
the solar spectrum for all elements except Li, S, K, Ir, Th, 
and U, so that the differential abundances are not strongly 
influenced by splitting effects. 

This is true if the lines have approximately the same 
strength in both stars. However, HD221170 and the Sun 
are different stars and it is quite difficult to find lines of 
similar intensity. For example, for Hf the equivalent widths 
of two of lines used are of the order of 2 and 3.5 mA in the 
Sun and near 45 and 15 mA in HD221170 . The first fine, 
A 3918.094 A, with an equivalent width of 2 and 45 mA in 
the spectra of the Sun and HD221170 , respectively, shows 
a relative abundance difference of hafnium -0.97 dex in the 
atmosphere of HD221170 with respect to the solar atmo- 
sphere. In the case of the second line, A 4093.155 A, the 
difference in the equivalent width is smaller (3.5 and 15 
mA) and the relative abundance is -1.49 dex. We can ex- 
pect that taking into account the splitting of the lines will 
decrease the overabundance of Hf with respect to iron in 
the atmosphere of HD221170 . The value obtained using 
the second line is therefore expected to be more reliable. 
But even this value shows a non-negligible deviation with 
respect to the solar r-process pattern, as will be shown in 
the next section of this paper. 

For the lines of several heavy elements (U, Th, 
Os, Ir), no counterparts in the solar spectrum ex- 
ist, and therefore the latest values of the oscillator 



Fig. 2. Same as Fig. 1 in the vicinity of the Pb I 4057.806 
A line. The three different synthetic spectra correspond to 
the optimum lead abundance and a deviation by ±0.5 dex 
from this value. Note that the unmarked strong line on 
the blueside of the Pb line corresponds to the ^^CH line. 
The separation of the ^^CH and Pb lines is smaller than 
the instrumental profile, so that the uncertainties in the 
determination of the wavelength and oscillator strength 
can strongly affect the Pb abundance determination. 
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Fig. 3. Same as Fig. 1 in the vicinity of the Th II 4019.129 
A line. The different synthetic spectra correspond to 
the mean thorium abundance of -1.18 (in the scale log 
N(H)=12) and a deviation by ±0.5 dex from this value. 



strength (Nilsson et al. 2002a, b, Ivarsson 2003) are 
adopted. Several lines of hafnium, lead, thorium and 
uranium can be found in Fig. 1-6. The line data of 
iron and other chemical elements observed in the pho- 
tosphere of HD221170 can be obtained in table for- 
mat at the websites: "users . Odessa. net/""yua" and 
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Table 3. The mean abundances of chemical elements in the atmosphere of HD221170 with respect to their abundances 
in the solar atmosphere. 
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Fig. 4. Same as Fig. 3 in the vicinity of the Th II 4094.747 
A line. 



Fig. 5. Same as Fig. 3 in the vicinity of the Th II 5989.045 
A line. 



"yushchenko .netf irms . com" as well in electronic Tables 
1,2 at the website of this Journal. 

Table 1 and 2 contain data on lines for iron and 
other chemical elements. Table 1 provides the ioniza- 
tion stage, equivalent width, oscillator strength, energy of 
lower level and derived iron abundance for each iron wave- 
length considered. For the other elements. Table 2 gives 
in addition to the line data (identification, wavelength, 
oscillator strength, the energy of lower level), the abun- 
dance determination (relative abundance with respect to 
the solar system value, and absolute values of the abun- 
dance calculated from this line using the spectrum of 
HD221170 and the solar spectrum), the levels of blend- 
ing of the line in the synthetic spectra of HD221170 and 
of the Sun, the depths of the line in the synthetic spectra 
of HD221170 and of the Sun, and the abundances in the 
atmosphere of HD221170 , calculated with two different 
atmosphere models, namely for a surface gravity increased 
by 0.2 dex, and for an effective temperature increased by 
100 K. 

In Table 3, the mean elemental abundances in the at- 
mosphere of HD221170 are given with respect to the solar 
atmosphere value. For all investigated elements. Table 3 
includes information on the atomic number, designation 
of the ionization stage, as well as the relative abundances 
(the last digits in brackets correspond to the estimated 
errors) and the number of lines determined in Paper I 
and in the present study. Mean abundances are calculated 
with the best set of atmosphere parameters, with a surface 
gravity increased by 0.2 dex and with an effective temper- 
ature increased by 100 K. The absolute abundances are 
available in the electronic version of the table. A com- 
parison of our abundances with previously published data 
can be found in Paper I. The abundances of S and Pb 
can be fomid here for the first time. The abundances of 
Mo, Er, Tm, Hf, W, Os, Ir, Th, and U first determined 
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Fig. 6. Same as Fig. 1 in the vicinity of the U II 4050.041 
A line. The different synthetic spectra correspond to a 
uranium abundance of -2.52, -2.02, and -1.72 (in the scale 
log N(H)=12). 

by Yushchenko et al. (2002) are updated in the present 
study. 

In Table 4 we show the differential and absolute abun- 
dances of iron and some of the key elements used to 
trace back the nucleosynthesis and corresponding cos- 
mochronometry discussed below. The errors correspond 
to the standard deviations of abundances derived from the 
individual lines of the element. For uranium and thorium, 
only lines with new oscillator strengths are considered. 

Concerning the detection of carbon and nitrogen, their 
atomic lines are hardly detectable, so that their abundance 
determinations are usually based on molecular lines. We 
calculated several synthetic spectra in the whole observed 
region using different C and O abundances and found 
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Table 4. Abundance of iron and of some elements de- 
termined in the present work. The columns provide, re- 
spectively, the atomic number, the clement symbol, the 
number of lines analyzed, the mean abundance in the at- 
mosphere of HD221170 with respect to solar photosphere 
(Grevesse & Sauval 1998) and the abundance in the scale 
logN(H)=12. 



z 


Ident. 


N 


Alog N 


logN 


26 


Fe I 


221 


-2.09-f0.07 


5.41±0.07 




Fe II 


30 


-2.09+0.17 


5.41±0.17 


63 


Eu II 


6 


-1.54-F0.13 


-1.03±0.13 


76 


Os I 


4 


-1.23-F0.12 


0.22±0.12 


77 


Ir I 


1 




0.17±0.20 


82 


Pb I 


1 


-1.85-1-0.20 


0.10±0.20 


90 


Th II 


5 




-I.18i0.ll 


92 


U II 


5 




<-2.02±0.20 



■1.0 



o 



-1.5 



-2.0 



-2.5 



-3.0 




10 20 30 40 50 60 70 80 90 100 

Z 

Fig. 7. The abundances of chemical elements and ions 
in the atmosphere of HD221170 with respect to the solar 
atmosphere value. Triangles correspond to upper limits for 
Li, Al, U. 



abundances close to -2.7 and -2.13 with respect to the 
Sun. This result is in good agreement with Sneden et al. 
(1986) and Craft et al. (1992), so that it is adopted in the 
abundance determination of the other elements. Special 
attention was also paid to the determination of the iso- 
topic carbon ratio which can strongly influence the final 
results. To estimate the ^^C/^^C ratio, several synthetic 
spectra were calculated and compared with observation. 
We finally obtained the value of ^^C/^'^C=6 in agreement 
with the determination of Sneden et al. (1986). 

The wavelengths of several ^^CH lines were changed ac- 
cording to Johnson & Bolte (2001). The corresponding list 
of spectral lines in the vicinity of the Th II 4019. 128 A line 
was considered in our calculations. Note that our result 
concerning the thorium abundance is based on 5 lines with 
Nilsson et al. (2002) oscillator strengths. Similar abun- 
dances are obtained with all these lines. We therefore 



recommend them for future abundance determinations in 
other stars. 

4. Discussion 

About half of the stable nuclei heavier than iron in the 
Universe are synthesized by the rapid neutron-capture 
process, also known as the r-process. The r-process is 
believed to occur predominantly in the latest stages of 
evolution of massive stars (heavier than about 10 solar 
masses) during their supernova explosion. Fig. 8 compares 
HD221170 abundances with the solar content in r-process 
nuclei. It clearly shows that for elements heavier than Ba 
(Z=56), the HD221170 surface abundance pattern is very 
similar to the solar one. Globally, it confirms the previous 
observation of ultra-metal-poor stars (Westin et al. 2000; 
Cayrel et al. 2001; Hill et al. 2002; Cowan et al. 2002; 
Sneden et al. 1998, 2003) that already in the early age of 
the Galaxy, the r-process was operational and quite unique 
in its production of, at least, the Z=56 to 78 elements all 
along the life of the Galaxy. However, the universality of 
the r-process for the production of elements Z>56, includ- 
ing the Pb-peak elements and the actinides, still remains 
to be confirmed. Such a universality has deep implica- 
tions including in particular the invariance of the relative 
r-nuclidic abundances with respect to galactic chemical 
evolution effects and the possibility to develop a stellar 
chronometry based on the actinide content of the metal- 
poor stars and consequently to estimate a lower limit to 
the age of the Galaxy. 

Although for most of the elements with Z=56-78 the 
present observations seem to comfirm this conclusion, the 
specific cases of Eu and Hf observed in HD221170 show 
the first indication of a non-negligible deviation with re- 
spect to the solar pattern (Fig. 8). Such a deviation re- 
mains rather small, but, as shown by Goricly & Arnould 
(1997), large deviations are not expected in this mass re- 
gion because of the nuclear correlations inherent to the 
nuclear aspects of the r-process nucleosynthesis. It was 
shown in the previous section that taking into account the 
hyperfine and isotopic structure of Hf lines could reduce 
the deviation but not to the extent to become compati- 
ble with the solar pattern. In the case of Eu, the abun- 
dance determination includes the splitting effects in both 
HD221170 and the Sun. 

In this comparison with the solar r-abundance distri- 
bution, another interesting feature of HD221170 is its high 
Pb surface abundance. Such an unusually high abundance 
has been already observed in the previous ultra-metal- 
poor r-process-enriched stars except CS31082-001. Lead 
can be efficiently produced by the r-process, but also by 
the slow neutron-capture process, known as the s-process, 
during the AGB phase of low-mctallicity stars with masses 
in the range 0.8-8 solar masses (Van Eck et al. 2001). Such 
a Pb s-enrichment is however systematically accompanied 
by a simultaneous carbon enrichment. For this reason, the 
significant Pb enrichment of a star like CS22892-052 char- 
acterized by a large C/0>1 ratio could be explained by 
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Table 5. Comparison of abundance ratios (in log scale) 
in Th-rich stars. The second column gives the metallicity 
in terms of [Fe/H]. 



Star 


[Fe/H] 


Th/Eu 


Th/Pb 


Th/U 


HD221170 


-2.1 


-0.15 


-1.28 


>0.84 


0831082-001" 


-2.9 


-0.22 


>-0.78 


0.94 


CS22892-052'' 


-3.1 


-0.62 


-1.62 


>0.73 


HD115444'= 


-3.0 


-0.60 


-1.8"* 


>0.40 


BD-M7°3248'= 


-2.1 


-0.51 


>-1.48 


>0.82 



"Hill et al. (2002). 
''Sneden et al. (2003). 

''Wostin ct al. (2000) for all ratios except Th/Pb. 
■^Th: Westin et al. (2000), Pb: Sneden et al. (1998). 
•^Cowan et al. (2002). 



additional C and Pb s-process enrichment from a com- 
panion star. But, in the specific case of HD221170 , its 
low C/O=0.025 ratio characteristic of heavy mass stars 
tend to argue against this possibility. The high Pb con- 
tent of HD221170 provides the second strong hint that 
the r-proc;c!ss might not be universal in its production of 
the elements above Ba. 

Finally, the present observation provides an accurate 
determination of the surface Th as well as an upper limit 
of U abundances in HD221170 . These actinides have been 
used extensively in recent years to estimate the age of the 
oldest metal-poor stars of the Galaxy and in so doing, a 
lower limit to the age of our Galaxy. The Th/Eu, Th/Pb 
and Th/U abundance ratios in HD221170 are compared 
in Table 5 to the corresponding ratios determined in the 
other four r-process-rich metal-poor stars. These ratios are 
of particular interest since they compare the Th abun- 
dance to light species (represented here by Eu), Pb (be- 
lieved to be produced in the same "environment" as Th, 
see e.g Goriely & Clerbaux, 1999) and U, a neighboring 
element. Based on the universality of the r-process, a rel- 
atively reliable Th/U chronometry can be built. An abso- 
lute determination of the Th and U abundances has been 
achieved for CS31082-001 only, leading to an age estimate 
of 13±4 Gyr (Hill et al. 2002; Goriely & Arnould 2001). 
Prom Table 5, it can be deduced from the Th/U ratio that, 
within the universality assumption, HD221170 is older or 
at most 2 Gyr younger than CS31082-001. 

The HD221170 Th/Eu ratio is seen to be the highest 
among the different stars, which is compatible with its rel- 
atively higher metallicity [Fe/H]. However, assuming the 
universality of the r-process, i.e the same original actinide 
production in all these stars (relative to the Z = 56 — 78 
r-elements), this would imply that HD22117G is younger 
than CS22892-052 by about 22±4 Gyr, but also younger 
by 17±4 Gyr than BD-hl7°3248, a star of similar metaUic- 
ity. This would imply that BD-|-17°3248 would be at least 
28 Gyr old, an observation that is difficult to reconcile 
with the traditional Big-Bang cosmology or the WMAP 
observation It is also seen that as far as the Th/Eu ra- 
tio is concerned, HD221170 is compatible with CS31082- 




40 50 60 70 80 90 

z 

Fig. 8. Comparison of the surface abundances in 
HD221170 with the solar system r-abundance distribution 
scaled to the observed Er abundance. The solar system 
abundances corresponds to the photospheric determina- 
tion of Grevesse & Sauval (1998) with the solar r-process 
distribution estimated by Goriely (1999). 

001, clearly confirming that the Th-rich star CS31082-001 
is not a rare exceptional case. This conclusion was also 
drawn by Honda et al. (2004) who determined at the sur- 
face of CS30306-132 a ratio of [Th/Eu]=-0.10 even higher 
than for CS31082-001 or HD221170. The Th/Pb ratios 
shown in Table 5 again confirms the previous conclusion 
that there is a large scatter in the r-abundances for the 
heaviest elements. The Th/Pb chronometry assuming the 
universality of the r-process would lead to the conclu- 
sion that HD115444 is older than CS31082-001 (of similar 
metallicity) by about 48±14 Gyr. 

5. Conclusion 

We determined the abundances of 43 chemical elements 
(including Th and U) in the atmosphere of the bright halo 
star HD221170. The present spectroscopic observation in- 
dicates that the r-process is not universal or more exactly 
the astrophysical site(s) hosting the r-process do(es) not 
always lead to a unique relative abundance distribution for 
the bulk Ba to Hg elements, the Pb-peak elements, and the 
actinides. Although the previous observations of CS22892- 
052, HD115444 and BD-M7°3248 suggested the universal 
feature of the r-process, the present HD221170 analysis 
confirms the non-universal trend originally suggested by 
the CS31082-001 observation. This new finding rules out 
the use of thorium chronometry to estimate the age of the 
oldest stars in the Galaxy. 
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